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It is theoretically shown that storage and time-reversal of light pulses can be achieved in a coupled- 
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resonator optical waveguide by dynamic tuning of the cavity resonances without maintaining the 
translational invariance of the system. The control exploits the Bloch oscillation motion of a light 
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pulse in presence of a refractive index ramp, and it is therefore rather different from the mechanism 

O ■ 

of adiabatic band compression and reversal proposed by Yanik and Fan in recent works [M.F. Yanik 
•Sh . and S. Fan, Phys. Rev. Lett. 92, 083901 (2004); Phys. Rev. Lett. 93, 173903 (2004)]. 
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The possibility of dynamically control the resonant properties of microresonator systems 
via small refractive index modulation represents a promising and powerful approach for an 
all-optical coherent control of light in nanophotonic structures jl, 2,1^,0, [5]. Recently, several 
theoretical papers have shown that a temporal modulation of the refractive index in pho- 
tonic crystals (PCs) and coupled-resonator optical waveguides (CROWs) can be exploited 
to coherently and reversibly control the spectrum of light, with important applications such 
as all-optical storage of light pulses , time reversal [2J and wavelength conversion 
The existence of a frequency shift on the spectrum of a light pulse reflected by a shock-wave 
front traveling in a PC was first pointed out by Reed et al. B], and adiabatic wave- 
length conversion by simple dynamic refractive index tuning of a high-Q microcavity in a 
PC has been numerically demonstrated in Ref.jlj]. In a series of recent papers, Yanik and 
Fan showed that an adiabatic and translationally-invariant tuning of the refractive index in 
a waveguide-resonator system can be exploited to stop, store and time-reverse light pulses 

ll, I2J, y, 19J . The general conditions requested to coherently stop or reverse light pulses have 
been stated in Refs. , and the possibility of overcoming the fundamental bandwidth- 
delay constraint of static resonator structures has been pointed out. The basic idea of these 
previous papers is that the band structure of a translational-invariant waveguide-resonator 
system can be dynamically modified by a proper tuning the refractive index without breaking 
the translational invariance of the system. For instance, stopping a light pulse corresponds 
to an adiabatic band compression process: an initial state of the system, having a relatively 
wide band to accommodate the incoming pulse, adiabatically evolves toward a final state 
in which the bandwidth shrinks to zero Q]. In practice, the adiabatic evolution is attained 
by a slow change of the refractive index of certain cavities forming the photonic structure 

l], [2], [J. The condition that the dynamic refractive index change does not break the trans- 
lational invariance of the system is important because it ensures that: (i) the system can 
be described in terms of a band diagram with a dispersion relation uo = u(k) relating the 
frequency u and the wave vector k of its eigenmodes; (ii) different wave vector components 
of the pulse are not mixed, so that all the coherent information encoded in the original pulse 
are maintained while its spectrum is adiabatically changed 4, 9]. 

In this work it is shown that a coherent and reversible control of light in a photonic structure 
by dynamic refractive index change does not necessarily require to maintain the translational 
invariance of the system. We illustrate this by demonstrating the possibility of stopping and 
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time-reversing light pulses in a CROW 
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12] with a dynamic refractive index gradient. 



In this system, light stopping and reversal is not due to adiabatic shrinking and reversal of 
the waveguide band structure, as in Refs.jl, 2], but it is a consequence of the coherent Bloch 



oscillation (BO) motion of t 
that, thought temporal [T" 



re light pulse induced by the index gradient. It is remarkable 
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15| and spatial [16l . Il71 ] BOs and related phenomena have 



been studied to a great extent in several linear optical systems, they have been not yet 
proposed as an all-optical means to stop or time-reverse light pulses. 

We consider a CROW made of a periodic array of identical coupled optical cavities, and 
indicate by u n = ujq + 5u n (t) the resonance frequency of the n-th cavity in the array, where 
8uj n {t) is a small frequency shift from the common frequency Uq which can be dynamically 
and externally changed by e.g. local refractive index control, as discussed in previous works 
[l, 5]. Practical implementations of CROW structures have been demonstrated in photonic 

jjl . ]j| or in a chain of coupled microrings [jjj]. In 



crystals with coupled defect cavities 
most cases, coupled mode theory [l|, 



11. 



201 ] can be used to describe the evolution of the 



field amplitudes a n in the cavities and therefore the process of coherent light control; the 
results obtained from coupled-mode theory have been shown in fact to be in excellent agree- 
ment with full numerical simulations using finite-difference time-domain methods (see, for 
instance, . For our system, coupled-mode equations read 



dt 



-«(a n _i + a n+ i) - 5u n {t)a r , 



(1) 



where k is the hopping amplitude between two adjacent cavities, which defines the bandwidth 
(4k) of the CROW. Note that cavity losses are not included in Eqs.(l), however a non- 
vanishing loss rate would just introduce a uniform exponential decay in time of a n which 
would set a maximum limit to the achievable delay time, as discussed in Ref.jl]. As in 
Refs. , field propagation is considered at a classical level; a full quantum treatment, 

which would require the introduction of noise sources in Eqs.(l) to account for quantum 
noise, is not necessary for the present analysis which deals with passive CROW structures. 
Contrary to Refs. []], we assume that the modulation of cavity resonances used for 
coherent light control is not translational invariant, i.e. 5u n depends on n. Precisely, we 
assume that a ramp with a time- varying slope a(t) is imposed to the resonances of N 
adjacent cavities in the CROW, leading to a site-dependent frequency shift 5uj n {t) = na(t) 
for 1 < n < N and Soo n (t) = for n > N and for n < 1. Note that the dynamic part 
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of the CROW which realizes stopping or time reversal of light is confined in the region 
1 < n < N, which is indicated by a rectangular dotted box in Fig.l. The total length of 
the system realizing stopping or time reversal of light is therefore L = Nd, where d is the 
distance between two adjacent cavities. The modulation a(t) is assumed to vanish for t < t\ 
and t > t 2 [see Fig. 1(a)], [ti,t 2 ] being the time interval needed to stop or time-reverse an 
incoming pulse. The switch-on time ti is chosen just after the pulse, propagating along the 
CROW and coming from n = — oo, is fully entered in the dynamic part of the CROW, 
whereas the length L is chosen long enough to ensure that the pulse remains fully confined 
in the cavities 1 < n < N for the whole time interval [ti,^]- Note that, as for t < t\ and 
t > t 2 the pulse propagates in the CROW at a constant group velocity, during the time 
interval t 2 ] the pulse motion as ruled by Eqs.(l) is more involved and turns out to be 
fully analogous to the motion of a Bloch particle, within a tight-binding model, subjected to 



a time-dependent field a(t) (see, e.g., [2l|, |22j, |23j). It is indeed such a Bloch motion that can 



be properly exploited to stop or time- reverse an incoming pulse. In fact, let us suppose that 
the incoming light pulse, propagating in the forward direction of the waveguide and coming 
from n = — oo, has a carrier frequency tuned at the middle of the CROW transmission band 
and its spectral extension is smaller than the CROW band width 4k. For t < t\ one can 
then write (see the Appendix for technical details) 

dQ F(Q)exp(iQn + 2iKt cos Q), (2) 

-7T 

where the spectrum F(Q) is nonvanishing in a small region at around Q = Qo = 7r/2. 
The shape of F(Q) is determined from the excitation condition of the CROW at n — > — oo 
or, equivalently, from the field distribution a n (io) along the CROW at a given initial time 
t = t < t\. For instance, assuming without loss of generality t = 0, the latter condition 
yields for the spectrum F(Q) = l/{2iz) J2 n fl n(0) exp(— iQn) (see the Appendix). For to < 
t < ti, Eq.(2) describes a pulse which propagates along the waveguide with a group velocity 
v g = 2dKsinQ Q = 2dn. At time t — ti, we assume that the pulse is fully entered in the box 
system of Fig.l, and the modulation of cavity resonances is then switched on. The exact 
solution to Eqs.(l), which is the continuation of Eq.(2) at times t > t\, can be calculated in 
a closed form and reads (see the Appendix) 

dQ F(Q)exp[iQn + i9(Q,t)}, (3) 

-7T 
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FIG. 1: Schematic of a dynamic CROW with a linear gradient of resonances, (a) Distribution of 
cavity resonances (left) and temporal behavior of gradient amplitude a(t) (right); (b) Process of 
pulse storage; (c) Process of time reversal. 

where we have set ^(t) = f* dt'a(f) and 9(Q, t) = 2k f* dt' cos[Q + 7(t')] + 2kI\ cos(Q). At 
time t 2 = t\ + t, the modulation is switched off, and for t > t 2 one then has 



a n (t) = exp(i7 n) / dQ F(Q) exp[if(Q) + i<f>(Q)} x 

J —TT 

x exp[iQn + 2intcos(Q + jo)], (4) 



where we have set 7 = f t ^ dt a(t), f(Q) = 2/ttxcosQ — 2nt 2 cos(Q + j ), and 
4>{Q) = QiQiti) = 2« Jj* 2 dt cos[Q + j(t)]. The modulation parameters are chosen to either 
store or time reverse the incoming pulse (see Fig.l). In both cases, we assume that the 
length L of the system is large enough to entirely contain the pulse in the whole interval 
[ti,t 2 ]- in case of pulse storage, after the modulation is switched off the pulse escapes 
from the system in the forward direction with the same group velocity v g = 2dn as that 
of the incoming pulse, but it is delayed by a time ~ r [see Fig. 1(b)]. In case of time 
reversal [Fig. 1(c)], the incoming pulse is reflected from the system, which thus acts as a 
phase-conjugation mirror. 



Consider first the process of pulse storage. To this aim, let us assume that the area 
7o be an integer multiple of 2ir. In this case, for t > t 2 from Eq.(4) one obtains 

/TT 
dQ F(Q) exp[i(j)(Q)] exp[iQn + 2in(t — r) cosQ]. (5) 
■TT 

A comparison of Eqs.(2) and (5) clearly shows that, if <p = the effect of the modulation is 
that of storing the pulse for a time r = t 2 — t\ without introducing any distortion: in fact, 
one has a n {t2) = a n (t\). If the area in an integer multiple of 2u but 7^ 0, the additional 



5 




(b) 



(c) 



20 40 60 80 
Propagation distance zld 



100 °o 



0.05 0.1 
Gradient a(f)/K 



c 0.8 

■S 0.6 

£ 0.4 
'co 

§ 0.2 



2 3 4 

Normalized Time ?/?„„ 



FIG. 2: (color online) Storage of a Gaussian pulse in a coupled resonator waveguide consisting of 
N = 100 cavities, (a) Gray-scale plot showing the space-time pulse intensity evolution (note the 
characteristic BO motion), (b) Profile of the applied modulation a(t) = an exp{— [(nt — 150)/to] 6 } 
with ao/n = 0.0958 and tq = 106, corresponding to an area 70 = Qtt. (c) Process of pulse storage: 
the solid black curve is the intensity profile of the incoming pulse as recorded in the first cavity 
of the waveguide (z = 0), whereas the solid red curve and dashed black curve are the intensity 
profiles of the outcoming pulse as recorded in the last cavity of the waveguide (z = L = lOOd) in 
the presence and in the absence of the modulation, respectively. In (c) time is normalized to the 
transit time t pass = N/(2n) = 50/k of the pulse in the system. 

phase 4>(Q) may introduce a non-negligible pulse distortion. The distortionless condition 
cf) = is exactly satisfied in two important cases: a step- wise modulation a(t) = an const 
(with a r = 2nl, I is an integer), and a sinusoidal modulation a(t) = ancos(fit), with 
tQ = 27rl and Jo(ao/Q) = 0. These two cases realize the well-known dc or ac BO motion 
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23| of the light pulse in the interval [ti , t 2 ] : pulse storage is therefore due to the 



periodic motion of the light pulse which returns to its initial position after each BO (or ac 
field) period. It is worth commenting more deeply the very different mechanisms underlying 
light stopping in the translational- invariant system of Ref . [l| with the one considered in the 
present work. In Ref. the modulation of cavity resonances preserves the translational 
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FIG. 3: Time reversal of a non-symmetric double-peaked optical pulse in a coupled resonator 
optical waveguide consisting of N = 60 cavities, (a) Gray-scale plot showing the space-time pulse 
intensity evolution, (b) Profile of the applied modulation a(t) = anexp{— [( K t — 50)/rn] 6 } with 
oto/n = 0.7 and tq = 17, corresponding to an area 70 = 7ir. (c) Pulse intensity profile versus scaled 
time nt as recorded in the first cavity of the waveguide (z = 0). 



symmetry and, as a consequence, cross talk between different wave vector components of 
the pulse is prevented as the waveguide band accomodating the pulse shrinks to zero and 
the pulse group velocity adiabatically decreases. Additionally, the tuning process must be 
slow enough to ensure reversibility, i.e. to ensure the validity of the adiabatic theorem. In 
the present work, the tuning of the cavity resonances breaks the translational symmetry 
of the system and the different wave vector components Q of the pulse undergo a drift 
motion in the reciprocal space according to the well-known 'acceleration theorem' of a 
Bloch particle studied in solid-state physics [3] (see the Appendix for more details). The 
motion of Q in the reciprocal space is accompanied by a shift of the pulse carrier frequency, 
which spans in a periodic fashion the full band of the waveguide, and by a periodic motion 
of the pulse in the "stopping box" of Fig.l, which is hence trapped inside it [see Fig. 2(a) 
to be discussed later]. In particular, for a step-wise modulation a(t) which will be mainly 
considered in this work 24J, the temporal periodicity of the motion is tb = 2n/ao. This 
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is the well-known periodic Bloch motion which is related to the existence for Eqs.(l) of a 
discrete Wannier-Stark ladder spectrum instead of a continuous band spectrum (for more 



22] ). Therefore, as in the translational- invariant waveguide system 



details see, for instance, 

of Ref. [l| light stopping is achieved by adiabatically shrinking to zero the band of the pulse, 
in our system light stopping can be viewed as a trapping effect due to the appearance of 
the periodic Bloch motion in the dynamic part of the CROW structure. Note that, as 
opposed to the method of Ref. in our case adiabaticity of the tuning process is not 
required, however the stopping time r is quantized since it must be an integer multiple of 
the Bloch period tb- Nevertheless, with a suitable choice of the gradient «o (and hence 
of tb), a target delay time r can be achieved. If t p is the duration of the incoming pulse 
to be delayed (with t p < we can estimated the minimum length L of the system as 

L = L p + Lb, where L p ~ t p v g = 2dnT p is the spatial extension of the pulse in the waveguide 
in the absence of the modulation and ~ And/ao is the amplitude of the BO motion. 
Hence the minimum number of cavities of the system is given by N = L/d ~ 2k[t p + 2/ao). 
It should be noted that in practice a sharp step-wise modulation can never be realized, 
and a finite rise time during switch on and off should be accounted for. Though <p(Q) does 
not exactly vanish in this case, dispersive effects can be kept however at a small level. As 
an example, Fig.2 shows the process of light storage as obtained by a direct numerical 
simulation of Eqs.(l) using a super-Gaussian profile for the gradient a(t). The area 70 is 
chosen to be 6n, so that pulse trapping corresponds to three BO periods, as clearly shown 
in the space-time plot of Fig.2(a). Note that the system comprises N = 100 cavities, and 
therefore the length L of the waveguide needed to perform light storage is L = Nd = lOOd. 
The intensity profile of the nearly Gaussian-shaped incoming pulse in the initial cavity 
(z = 0) is indicated by the black solid line in Fig.2(c). In the figure, the intensity profile 
of the outcoming pulse at the last cavity of the system (z = 100c?) is shown by the red 
solid line, whereas the dashed curve indicates the intensity profile of the output pulse at 
the last cavity in the absence of index gradient, i.e. when the pulse freely propagates along 
the system at the group velocity v g = 2nd. In Figs. 2 (a) and 2(b), time is normalized to 
l/n, whereas in Fig.2(c) time is normalized to the transit time t pass = L/v g = N/(2k) 
of the pulse in the system. Note that the maximum frequency shift of cavity resonance 
needed to achieve the process of pulse storage is Su max = ±(N/2)ao ~ ±5k. Assuming 
that a change 5n of the refractive index n produces a change 5uj ~ Uo(5n/n) of the cavity 



S 



resonance uq, the index ramp of Fig. 2(b) thus corresponds to a maximum refractive index 



ight stopping 
1]. To get an 



change 5n/n ~ 5k/uq. This value is comparable to the one requested for 
by means of adiabatic band compression in the translational-invariant case 
idea of typical values in real physical units, let us assume e.g. a carrier angular frequency 
Uq ~ 1.216 x 10 15 rad/s (corresponding to a wavelength A ~ 1.55 /im) and a maximum index 
change Sn/n ~ 5 x 10~ 4 , which is comparable to the one used in previous studies (see, e.g. 



, |5j). The bandwidth 2k of the waveguide and the transit time t pass in the figure are then 
given by 2k ~ 2 x KTV, ~ 2.43 x 10 11 rad/s (i.e. ~ 39 GHz) and t pass = N/(2k) ~ 410 
ps, respectively. For such parameter values, Fig.l simulates the stopping of a ~ 68 ps-long 
(FWHM) Gaussian pulse with a storage time r ~ 1.75 ns. 



The process of time-reversal of a light pulse is simply achieved when the area 70 is 
equal to 7r, apart from integer multiples of 2tt. In fact, in this case for t > t 2 from Eq.(4) 
one obtains 

a n (t) = (-l) n dQF(Q)exp[i(j)(Q)] x 

J — n 

exp[iQn — 2in{t — t\ — t 2 ) cosQ], (6) 

A comparison of Eqs.(2) and (6) clearly shows the sign reversal of the frequency 2k cos Q for 
any wave number Q in the integral term, which is the signature of time reversal of the pulse. 
Physically, time reversal is due to the fact that, for the it area, the spectrum of the wave 
packet in the reciprocal Q space (quasi-momentum) is shifted in the Brillouin zone from 
Q = n / 2 to Qq = — 7r/2, thus producing spectral inversion. In addition, since the group 
velocity is correspondingly reversed, the pulse is reflected by the system and thus propagates 
backward. As in the previous case, the process of time reversal does not introduce pulse 
distortion provided that = 0. For a step-wise modulation, contrary to the 2tt area case 
4>{Q) does not vanish and one has <f>{Q) = — (4/c/ao) sinQ. However, this additional phase 
term may be kept small by choosing e.g. a sufficiently large value of do, thus minimizing pulse 
distortion. An example of time reversal of an asymmetric pulse with minimal distortion, as 
obtained by a direct numerical simulation of Eqs.(l) using a super-Gaussian profile for the 
gradient oe(t), is shown in Fig. 3. Note that in this case the pulse undergoes a semi-integer 
number of BO periods. 

In conclusion, it has been theoretically shown that storage and time-reversal of light can be 
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realized by exploiting BOs in a dynamic coupled-resonator waveguide. The proposed scheme 
is rather distinct from the adiabatic band compression technique recently proposed in Refs. 

, and provides a noteworthy example of coherent light control in a system with broken 
translational invariance. 



APPENDIX A 

In this Appendix we provide a detailed derivation of the solution to the coupled-mode 
equations (1) at times t < t±, t\ < t < i 2 an d t > t 2 presented in the text ]Eqs.(2), (3) 



and (4)]. To this aim, we follow a rather standard technique (see, for instance, 21J) and 
introduce the time-varying Fourier spectrum G(Q,t) defined by the relation 

G(Q,t) = — Y] a n (t)exp(-iQn). (Al) 

Z7T z — ' 

n=— oo 

The amplitudes a n (t) can be derived from the spectrum G(Q,t) after inversion according to 
the relation 

dQG(Q,t)exp(iQn). (A2) 

Using the coupled-mode equations (1) with 6u n (t) = a(t)n, the following differential equa- 
tion for the spectrum G can be easily derived [25 ] 



dG OG „ , . — y . — ^ , . 

— + a— = 2ikGcosQ. A3 
ot oQ 

For t < ti, we have a(t) = 0, and therefore the solution to Eq.(A3) is simply given by 

G(Q, t) = F(Q)exp(2iKtcosQ) (t < h), (A4) 

where the profile F(Q) is determined by the spectrum G at an initial time t = to by means 
of Eq.(Al) once the field distribution a n (to) is assigned. Note that substitution of Eq.(A4) 
into Eq.(A2) yields Eq.(2) given in the text. 

For t > ti, a(t) is nonvanishing and the solution to Eq.(A3), which is a continuation of 
Eq.(A4) for times t > t%, can be easily obtained after the change of variables rj = t and 
£ = Q — 7(i), where we have set 

1 (t) = f dt'ait'). (A5) 
Ju 
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With these new variables, Eq.(A3) is transformed into the equation 

= 2ikG{^ n) cos[£ + 7(7/)], (A6) 

which can be easily integrated with the initial condition G(£, r\ = t\) = F{^) exp(2mii cos£). 
Upon re-introducing the old variables Q and t, one then obtains 

G(Q,t) = F(Q-7(t))exp[2i/cticos(Q-7(t))]exp [ 2iK J dt ' C0S lQ + 7(*') ~ 7(*)] } (* > *0 

(A7) 

Substituting Eq.(A7) into Eq.(A2), after the change of integration variable Q' = Q—j(t) and 
taking into account the 27r-periodicity of the spectrum G(Q,t) with respect to the variable 
Q, one then readily obtains Eq.(3) given in text. Note that \G(Q,t)\ 2 = \F(Q — 7(t))| 2 , i.e. 
the role of the gradient a(t) is to induce a rigid drift of the initial spectrum, a result which 



is known as " acceleration theorem" in the solid-state physics context 22J . In particular, for 
a constant gradient a(t) = ao, the drift of the spectrum is uniform in time. In this case, 
from Eq.(A5) it follows that, after a time r# = 2-n;/a$ from the initial time t = ti, one has 
r y(t 1 + T B ) = 2tt and, from Eq.(A7), G(Q,ti + r B ) = G(Q,ti), i.e. the initial field distribution 
in the CROW structure is retrieved: t b plays the role of the BO period which is determined 
by the gradient a . 

For t > t2, one has a(t) = and the spectrum G(Q,t) is still given by Eq.(A7), where 
according to Eq.(A5) one has j(t) = f t * dt'a(t') = 70 for t > t2- Note that for t > t<i 
Eq.(A7) can be cast in the following form 

G{Q,t) =G(Q,t 2 )exp[2m(t-t 2 )cosQ] (t > h) (A8) 

Substitution of Eq.(A8) into Eq.(A2), with G(Q,t2) given by Eq.(A7) with t = t 2 , and after 
the change of variable Q' = Q — 70 in the integral of Eq.(A2), one finally obtains Eq.(4) 
given in the text for the solution at times t > t 2 . 
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down and stopping in this case is thus analogous to adiabatic band compression of Ref.[l[. 
The requested modulation frequency f2, however, turns out to be too high (larger than k) to 
be reasonably achieved in practice. 
[25] In deriving Eq.(A3) we assumed that the relation co n = na(t) holds even for n < 1 and n > N, 
though in practice the dynamic tuning of cavity resonances is applied solely for the cavities 
with index n in the range 1 < n < N (see Fig.l). However, taking into account that a(t) 
vanishes outside the interval [t\, £2] and that the size N of the stopping box is large enough to 
fully contain the pulse during the delay interval, i.e. a n (t) ~ for n < 1, n > N in the time 
interval t\ < t < t2, the previous assumption does not change the solution of the problem. 



13 



